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Abstract 



Charged black hole solutions with pion hair are discussed. These can be 
used to study monopole black hole catalysis of proton decay. There also exist 
multi-black hole skyrmion solutions with BPS monopole behaviour. 
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I. INTRODUCTION 



It is known that baryons can emerge as solitons in an effective meson field theory. In these 
models the baryon number is identified with a topological charge. Since this identification 
was originally made by T.H.R.Skyrme 40 years ago the soliton solution is called the 
skyrmion. 

Skyrmion models provide a useful way of analysing the monopole catalysis of proton 
decay, as originally pointed out by Callan and Witten [Q . The decay process can be allowed 
because, in the presence of a monopole, the baryon number is no longer equal to the topo- 
logical charge of the meson field. In fact, there exist non-topological solitons with non-zero 
baryon number which can decay without topological problems. 

Skyrmions have also been used to examine the absorbtion of baryons by microscopic 
black holes @-f|. Solutions representing skyrmions partially absorbed by Schwarzschild 
black holes provide a semiclassical framework to study the absorption rate of a proton by 
a black hole of comparable size. However, this process is rather insignificant because black 
holes of the size of a proton have large fluxes of Hawking radiation which swamp the proton 
decay. 

In this paper we analyse black hole solutions with skyrmion hair and magnetic charge, 
providing the semiclassical framework in which to study monopole black hole catalysis of 
proton decay. The flux of Hawking radiation from these holes is less than ordinary holes 
and even vanishes in extremal cases J7[. The black hole mass for abelian monopoles has the 
lower bound 

M = = pm pl « 2.54 x lO -7 ^ , (1) 

where e is a rationalised electric charge in gaussian units determined from the Dirac quan- 
tisation condition (ep = 1) as p ~ 11.7 and m v \ is the planck mass which is approximately 
2.1768 x 10~ 8 Kg. This kind of small black hole may be considered to have been created in 
the early universe and remain as a stable relic today. More complicated monopole solutions 
are also possible in Grand Unified Theories, depending on the details of the Higgs sector ||, 
but we will restrict attention to the simplest case. 

In the extreme case, M = pm v \, we find the remarkable situation that multiple black hole 
solutions are possible in which the gravitational, electromagnetic, and strong forces between 
the monopoles are all in balance. In this respect the solitons behave in an analagous way to 
BPS monopole solutions in the Yang-Mills-Higgs system |J. 

We find that the non-topological skyrmion solutions are stable within the confines of 
our model. This is, in part, due to the fact that we have not included the electron. The 
black hole cannot swallow the proton whole because this would tip it over the extremal 
limit. Light, charged particles are therefore needed to carry away the proton's charge when 
it decays. We discuss a way of including this effect in the conclusion. 

II. THE LAGRANGIAN 

We shall consider models with a charged SU(2) meson field U. The Lagrangian is based 
on a gauged version of the original skyrmion Lagrangian || . It should be regarded as a min- 
imal form of the meson effective action, since extra terms could also be included. However, 
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some of the terms in the Lagrangian which couple the chiral field to electromagnetism can 
be deduced from current algebra techniques and were constructed by Callan and Witten || . 
The Lagrangian can be divided into into four parts, 

C = C x + £ 2 + £ 3 + U (2) 

where the meson parts are 

Ci = ||tr (u^D^UU^D^u) + ^tr ([U^D^U, U^D^U}) 

C 2 = ^e^ pa A^i{Q (dvUU^dpUU^dvUU- 1 + U^d.UU^dpUU^d^u)} 

?'e 2 / 1 1 \ 

+ ^e^{d,A u )A p ti \Q 2 d a UU~ l + Q^d.U + -Qd^UQU' 1 - -QUQd^U-'j 

and the free actions are 

£3 = -iib F - F " ^sb* < 3) 



The quantity Q is the charge matrix of quarks and D is a covariant derivative defined by 

2 
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Q= (4) 
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D fl U = d fi U-ieA fl [Q,U\. (5) 

The abelian gauge field A^ and electric charge e are in unrationalised units, and H = c = 1. 
In the spherically symmetric case with a magnetic charge the gauge field has the form 

A =p(l- cos 9)d<p + $dt (6) 

where p is a magnetic charge and the Dirac quantisation condition is pe = 1. The usual 
Skyrmion has a magnetic moment which would interact with a magnetic monopole and 
break the spherical symmetry. We use instead a non-topological ansatz for the chiral field, 

U = e if(r,t)* 3 _ ( 7 ) 

Despite the fact that this field is made up of neutral pions and commutes with the charge 
matrix Q, it has a non-zero total electric and baryonic charge due to the effects of anomalies, 
as we shall see in the next section. 



III. BARYON NUMBER AND ELECTRIC CHARGE 



The gauge invariant baryon current was constructed by Callan and Witten 

-— hru-^uu-^uu^dpu 

+3ieA l/ trQ(U- l d a UU~ 1 dpU - dJJU^dpUlJ- 1 ) 
+3ied 1/ A a trQ(U~ 1 d f5 U + dpUU' 1 ) 



Jb 



(8) 

(9) 
(10) 
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Only the third term survives for the chiral field anzatz (|7]), with a resulting baryon charge 



n B = f£ [f(oo) - f(0)} . (11) 

The solution with the boundary conditions /(0) = and f(oo) = 2tt posesses unit baryon 
number and hence it can be interpreted as a baryon surrounding the monopole. 

If the field anzatz (|7|) is substituted into the meson and electromagnetic interaction terms 
in the Lagrangian, they become 

A = ~%df)\ (12) 

C i = ~E i B i f, (13) 

C 3 = l(E>-B>), (14) 

where the index % = 1,2,3 and the electromagnetic fields E{ and Bi are defined by F$i = 
(—g tt ) l / 2 Ei and Fy = e ijk B k . When combined, 

The extrema of the action occur when the electric field is given by 

E = ^Bf (16) 

This situation is reminisent of the factorisation of the Lagrangian that occurs for a BPS 
monopole |J. 

The electric field implies a total charge 

or asymptotically 

goo = n B e (18) 

and the rig = 1 solution can therefore be interpreted as a proton. 

If a black hole appears in the background, the inner boundary condition for the field 
/ should be imposed not at the origin but at the event horizon r = r+. Thus the baryon 
number in the presence of an event horizon will be defined as 

n B = 2- [/(oo) - f(r + )\ . (19) 

If /(r+) = the baryon number is still an integer and conserved. This configuration repre- 
sents a proton tightly bound to the black hole. On the other hand if / takes some positive 
value at the horizon the baryon number is not an integer and the skyrmion carries fractional 
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baryon number and electric charge. This configuration will be interpreted as a proton par- 
tially swallowed by the black hole. In particular, f(r + ) = 2ir means that the black hole has 
swallowed a whole proton leaving nothing outside the horizon. 

It is interesting to observe that, while the baryon number disappears inside the horizon, 
the electric charge of the black hole can still be measured outside, turning the monopole 
black hole into a dyon black hole. Therefore, while the baryon number conservation is 
violated, charge conservation is not violated. 



IV. EXTREMAL BLACK HOLE SOLUTIONS 



In the extremal case we can obtain a general solution based on the Papapetrou-Majumbar 
metrics HT0| , p"Tf . We begin with the background metric fixed and later generalise to solve the 
full Einstein equations with chiral matter. The Papapetrou Majumbar metrics have the 
form 

ds 2 = -U' 2 dt 2 + U 2 (dx 2 + dy 2 + dz 2 ) (20) 

where 

U^ + t^ (2D 

n=l U n 

and R n is the ordinary Euclidean distance from the point mass M n located in three dimen- 
sional space. We also associate these point masses with magnetic charges P n = G x l 2 M n , 
and the magnetic field 

B = G- ll2 U~ x dU. (22) 
The matter lagrangian obtained earlier (|12"D has the form 

C = - 1 -F 2 (df) 2 -^B 2 (l + a 2 f 2 ) (23) 

where we will set 

a = f, /i 2 = nGF 2 (24) 

Z7T 

The skyrme field equation obtained from the Lagrangian on this background becomes 

-fi 2 d 2 f + a 2 U- 2 (dU) 2 f = (25) 

We loose no generality by taking equal charges P n = p. The solution with baryon number 
n B = n M is then 

/ = 2irU- s (26) 

where 
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(27) 

Since F w <C m ph we can use s ~ aju for the pion model. 

For a single black hole, the Reissner-Nordstrom coordinate r is related to R by R = r—r + 
and we obtain 

f = 2n(l- r -±) S . (28) 

The field is effectively expelled by the black hole and vanishes on the horizon r = r + . 

The mass of the chiral field configuration can be obtained by integrating the Lagrangian 



m f = l -Fl JjaifdS* (29) 
where £ is a large surface containing all of the masses. This gives 

m f = 2tt 3 sF 2 G J2 M n~ 7r 3/2 n B eF^ (30) 

n 

The total mass in the chiral field is much less than one baryon mass per mass point. We can 
see how it is energetically favourable for a free skyrmion to change its internal configuration 
from the original Skyrme form to the simpler form used here when it comes into contact 
with a black hole monopole. The topological description of this transformation for a single 
monopole is exactly as described in reference ||. 

It is interesting to see that the electrostatic energy cancells due to the factorisation 
occuring in the lagrangian flT5|). The chiral field mass is independent of the separation of 
the holes and therefore there are no forces between then. This is similar to situation for 
BPS monopole solutions ||, and suggests that there is a solution of the full Einstein-matter 
system. This existence of this solution will now be demonstrated. 

The spatial part of the Einstein tensor for the matric (p0[) is 



Gn = -ZU-^diU^U) + U- 2 (dU) 2 9ij (31) 

and the Ricci scalar is 

R = -2U~ l d 2 U (32) 



Substituting the Einstein tensor for the Lagrangian ( |23"D into the Einstein equations gives 

U~ l d 2 U = -fi 2 (df) 2 (33) 
U-^dJJXdjU) = -fi 2 W)W) + GB^il + a 2 f) (34) 

These make up a complete system of equations when we include the Maxwell equation 

d{UB) = (35) 

The second Einstein equation implies that df, dU and B are all parallel. We therefore 
impose a condition / = f(u), B = b(u)U~ 1 dU, where 
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u=-fi~ 1 logU (36) 

The system of equations becomes equivalent to an ordinary differential equation with inde- 
pendent parameter u, 

f" + /a (l + /«) f - a 2 b 2 f = (37) 

where 

1 4- f 12 

b2 = T^T2- (38) 

The horizon corresponds to u — ► — oo and the far region to u — > 0. The horizon must 
therefore be at a critical point of the first order system corresponding to (pTj). There is only 
one critical point, (/, /') = (0, 0), hence 

(/,/0-(0,0)as«--oo. (39) 

Since the critical point is a saddle, the solution is unique and exists for all values of \i. 
Having obtained the unique solution to (0), we then define 

V(u) = 1 + /X [° b{x)e~ tiX dx. (40) 

Ju 

It is easily seen from (|36|) that 

3 t V = V'diU = UB (41) 
Hence the Maxwell equation implies d 2 V = and we can write 

K = HE^ (42, 

Inverting ( fiO|) gives u(V). 

V. SPHERICALLY SYMMETRIC SOLUTIONS 

In the non-extremal case we shall consider spherically symmetric metrics which can be 
parameterised in the form 

ds 2 = -%e 2& dt 2 + r ^dr 2 + rW (43) 
r 2 A 

where A and 5 are functions of r and t. After inserting the metric and the other field ansatze 
(i)"(0) i n ^° the Einstein field equations, one obtains 

(AeV)' - = - ^-ti + £e-'f (44) 

5' = A (^- M / 2 + /°) (45) 

,46, 



where [l and A are constants, 

f = "fr. \' = (47) 

The electric charge within a sphere of radius r is given by equation (p!7|). 

For very small which is the case for pions, the chiral field has little effect on the 
background metric and we may take 5 = and express A in terms of the mass M, electric 
charge Q and magnetic charge p of the black hole as 

A = r 2 - 2GMr + G(Q 2 + p 2 ). (48) 

The skyrme field equation fl4"4|) on this background is therefore 

A 2 



(A/')' - 'jsf = 0, (49) 



This should be solved subject to the boundary condition on which fixes the total charge, 

Qoo = 7T"/oo = n B e. (50) 

The non-extremal black hole posseses two horizons at r = r_ and r = r + (r + > r_), related 
to the mass and charge by 

GM = ^(r_ + r + ), GQ 2 = r_r + - Gp 2 . (51) 

The solution to equation (EDT) can be obtained analytically, 



P (r±±r=) 

* \ r.\-—r- J 

where P s (z) is a Legendre function and 



/ = 2 ™b ^aX" (52) 



s = -- + 4/- + -^— (53) 
2 V4 r + r_ V ; 

The black hole becomes a dyon with electric charge related to the value of / at the event 
horizon, 

Q = - , r++ r-\ - ( 54 ) 



This relation can be solved, together with (|51|), to obtain Q = Q(M), showing the existence 
of a one parameter family of solutions (n# and p being regarded as fixed). In particular, 
Q — > as M approaches the extremal limit pm p i and the meson field is expelled from the 
hole. 
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Larger values of fi may be realised for a hypothetical model where U is unrelated to 
pions and this is discussed below. We will consider a static solution. We can replace A by 
a mass function m(r), defined implicitly by the relation 

A = r 2 - 2GW + G(p 2 + q 2 ) (55) 

where the charge is given by equation (|H]). The equations become 

w! = Aj' + ^ff (56) 
5' = A(/') 2 (57) 

Suitable boundary conditions are / — ► 2n and 5 — ► as r — > oo. 
In the numerical results the fields are scaled to the horizon size, 

t A Cm 
r = — , m = (59) 

The solutions are parameterised by a parameter p, defined by 

p =- r (60) 
r + 

which is restricted to p < 1. 

The extremal black hole solutions have A + = A' + = 0. The regular solution to equation 
(H) has, 

/+ = 0, Q = ^f + = 0, p = l. (61) 

Hence the proton lies fully outside the black hole, as we saw before. The numerical solution 
for / is shown in fig.(Q). This agrees well with the result P%| ) of the previuos section, because 
the value of [i used here is still quite small. The results are still qualitatively similar for 
chiral models with \x of order one. 
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FIGURES 
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FIG. 1. Skyrme field / as a function of f = r /r+ for an extremal hole and /j, = 10 4 . 

For the non-extremal solution, we begin the integration of the field equations close to 
the horizon, with 

in = rh + rhiif — 1) + m 2 (f — l) 2 + • • • (62) 
6 = 5 + + 5 1 (r-l) + --- (63) 
/ = / + + / 1 (f-l) + ... (64) 

where 5 + and /+ are shooting parameters determined so as to satisfy the boundary conditions 
at infinity. 

As can be seen from the above expansion, the skyrme field must have a nonzero value 
at the horizon otherwise the only allowed solution is the trivial one. Consequently the 
nonextremal black hole acquires an electric charge 

Q = £~f+, (65) 

and allows the skyrmion to have fractional electric charge. The numerical results for this 
solution are shown in fig.(|2|)-(|3|). Again, these agree well with the fixed background for small 
values of /x. 
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FIG. 2. Backreaction <5 (xlO 3 ) as a function of f = r/r + for a non-extremal black hole, 
p = pm p i/r + = 0.9. Results for fi = 1CP 3 and \i = 10~ 4 are shown. 
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FIG. 3. Mass function m as a function of f for p = 0.9. Results for [i = 10 -3 and \i = 10 -4 are 
shown 

We have a single one parameter family of solutions with p < 1. In fig.(^)-(^j), we plot the 
horizon radius r + and skyrmion mass mj as functions of black hole mass M. Figure (f|) is 
related to the entropy of the black hole (Anr 2 ^). The other figure shows how the proportion 
of the skyrmion which is swallowed by the black hole increases with the black hole mass. 



11 



r h /Mc 



4- 



1 1 1 1 1 1 M/Nfc 

0.0 0.5 1.0 1.5 2.0 2.5 3.0 



FIG. 4. Horizon radius r + /M c as a function of the black hole mass M/M c for \x = 10 
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FIG. 5. Skyrmion mass mj as a function of the black hole mass M/M c for \i = 10 4 . 

The last figure (§) shows how the horizon value of / changes as the coupling constant 
fj, changes. Since \x characterises the mass scale of the chiral model, this is amounts to a 
comparison of different models. As can be seen from the figure, for small \x the electromag- 
netic interaction is dominant so that the skyrmion is absorbed by the black hole to a lesser 
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extent. On the other hand for large \i the gravitational interaction is dominant so that most 
of the skyrmion is absorbed by the black hole. 




FIG. 6. The value of / at the horizon for various values of \x. 



VI. STABILITY ANALYSIS 



In this section we show that the skyrmion solutions which we have obtained are stable 
under spherically symmetric linear perturbations. We shall begin with the analysis of a 
skyrmion on the fixed background. 

In the fixed background case the skyrme field is the only perturbed field and can be 
expanded about the skyrmion solution fo by writing 



/M) = /oW+e'1(r). 
Equation flBTf) is inserted into equation to obtain the eigenvalue equation 



(66) 



A 



(67) 



where the background equation has been used. 

If to is real and u 2 > the solution is stable, and if to is imaginary and u 2 < it is 
unstable since the mode can grow or decay exponentially under the small perturbation. To 
show which is the case we multiply both sides of equation ( |67|) by £ and integrate in r from 
the horizon to infinity 



r+ 



dr 



to 



/*oo <r^ 
Jr+ IAq 



(68) 



13 



where integration by parts and boundary conditions have been used. It can be seen that the 
integrands of both sides are positive definite, which means that uo 2 > 0. Hence the skyrmion 
on the fixed background is linearly stable . 

Next we analyse the stability of the skyrmion with backreaction. In this case we have to 
expand the metric as well as the skyrme field around the classical solutions /n, 5q and A 

f(r,t) = f (r) + f 1 (r,t) 
6{r,t) = 5 + 5 1 {r,t) 
A(r,t) = Ao + Ai(r,t) . 



These are substituted into eq. (fP|)-(fP3D to obtain the following coupled equations up to first 
order 



/ A 2 



(AoJi/J + A f[ + Ai/jJ) e 5 ° - - (5 1 f + A) e 5 ° = -j-e^h 



An 



5[ = 2^rf Q f[ 
Ax 



Equation (ffl|) can be integrated with the help of the static field equation, 

A- L = -2AA0/0A • 



(69) 
(70) 
(71) 

(72) 



Substituting equation ( fT0|) and equation ([72]) into equation ( |6~9"D one obtains the first order 
equation for f\ 



(Aoe*/0'- 2 / . 2 (rA e 5 Vo 2 )' + 



A 2 (M ( 



fi = ^~ &0 h ■ 
^0 



(73) 



Setting fi(r,t) = £(r)e iu}t one obtains an eigenvalue equationfor £, 



(A e^')' + 2 / . 2 (rA e^ 2 )' + 



A 2 (M C ) 



^0 



(74) 



We introduce the tortoise coordinate r* such that 

dr* 1 



dr A e 5 ° 



(75) 



and r* runs from —00 to +00 as r runs from r + to +00. Then eq. (f74|) is reduced to the 
Sturm-Liouville equation 



dr* 2 



(76) 



where 



U 



2/. 2 (rAoeV^' + ^e* 



A e 5 °. 



(77) 
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On the left-hand side we have r 4 , which makes the equation different from the previous 
eigenvalue equation. However, since r 4 remains positive through the whole space, the same 
conditions for stability as the ordinary eigenvalue equation can be applied. As can be seen 
by examining U, U — > as r — > r + , (i.e. r* — > — oo), U — > > as r — > oo, and Z7 > 
in between. In addition the solution does not change its shape for any value of the coupling 
constant /i. Therefore we can safely conclude that a skyrmion with backreaction is also 
linearly stable. 



We are now able to describe some of the features of the interaction between a slowly 
moving proton and a black hole monopole. The free skyrmion has a magnetic moment and, 
if it has the correct orientation, it will be attracted to the monopole. When the proton 
approaches the black hole monopole, the fields rearrange themselves into the energetically 
prefered configuration of skyrmion hair solution described in this paper. 

The skyrmion solution on the fixed black hole background was obtained analytically in 
both extremal and nonextremal cases. In the extremal case, the baryon number and electric 
charge are expelled by the horizon and the system represents protons bound to monopole 
black holes. We found a general solution on the Eintein equations in the extremal case with 
many similarities to the BPS monopole system. 

In the nonextremal case, the monopole black hole partially swallows the proton and 
transforms to a dyon black hole, leaving fractional electric charge and baryon number outside 
the horizon. However, Hawking radiation cannot be ignored in this case. 

We have also obtained numerical skyrmion solutions in the nonextremal case with grav- 
itational backreaction. The effect of the backreaction is important only for very massive 
skyrmions i.e. very large coupling constant. The results are qualitatively similar to the 
solutions without backreaction. 

Since the solutions are stable, baryon decay can only take place when extra particle fields 
are included in the model. If we introduce a charged field on the fixed background, then 
following the proceedure described in reference ||, the equations become (approximately) 



The stability arguments no longer apply. The dynamical process of a black hole swallowing 
a proton can be examined by solving these time-dependent field equations numerically. 
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